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Abstract 

We present a computation of the inclusive Drell-Yan production cross-section in perturbative QCD to 
all orders in the limit of high partonic centre-of-mass energy. We compare our results to the fixed order 
NLO and NNLO results in MS scheme, and provide predictions at N^LO and beyond. Our expressions 
may be used to obtain fully resummed results for the inclusive cross-section. 
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1 Introduction 



Accurate perturbative calculations of benchmark inclusive cross-sections are an essential component 
of the LHC discovery programme. The most important of these benchmark processes are the Drell- 
Yan processes: production of jJ-'^-n" pairs, and the closely related vector boson production processes. 
Currently these cross-sections are known at NNLO in perturbative QCD [1-4]. The resummation of 
threshold logarithms is known up to N^LL [5-8]. 

When the invariant mass Q of the particles produced in the final state is well below the centre-of- 
mass energy, the typical values of x of the colliding partons may be rather small: x\X2 = Q^/S <^ 1. 
In fact this is true of most LHC processes: only when producing very massive states close to threshold 
do both partons carry a large fraction of the incoming longitudinal momentum. Whenever x is small, 
logarithms of x may spoil the perturbation series. Thus accurate calculations require the computation 
of the coefficients of these logarithms, and if they are large it may be necessary to resum them. 

The resummation of small- a; logarithms in the perturbative evolution of parton distribution functions 
at NLL is by now well understood (see for example rcf. [9] for a recent review). The general procedure 
for resumming hard cross-sections at the same order through /c^^-factorization is known [10, 11], and 
its implementation when the coupling runs understood [12]. Calculations have been performed for 
photoproduction processes [10,12], deep inelastic processes [11,13], hadroproduction of heavy quarks 
[10, 12, 14, 15], and gluonic Higgs production both in the pointlike limit [16], and for finite top mass 
[17, 18]. However for Drell-Yan and vector boson production the resummcd hard cross-section has yet 
to be computed. It is the purpose of this paper to perform this calculation, determining the coefficients 
of the leading high energy singularities of the Drell-Yan coefficient function to all orders in perturbation 
theory. This will enable the accurate evaluation of small-x corrections to these benchmark processes, 
soon to be measured at LHC. 



2 The Drell-Yan cross-section 

We wish to study the high energy behaviour of the Drell-Yan cross-section. We consider Uf quarks qi 
with electric charge Cj. The cross-section for the production of a lepton pair via an off-shell photon 
with squared momentum = can then be written as: 

<jirh,Q'') = MQ^) E T— ['—Dabi^,;asm)Fa{x,,Q'')F,{x2,Q^), (1) 

, _ Jp Xi J p X2 
a,b=qi,qj,g f f 

where = Q'^/S, and Fa{x, Q"^) = xfa{x, Q^) is the integrated parton density for parton a, the indices 
a and b running over the different initial partons qi, qj, g. The LO partonic cross-section is simply 



MQ') = ^2^y), (2) 



Q2 3iV< 



where (e^) = ^ ef : the dimensionless coefficient functions Dah then contain all the QCD radiative 
corrections. 
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(a) (b) (c) 

Figure 1: (a) LO and (b) NLO diagrams contributing to Dqq, and (c) a NNLO diagram contributing 
to the flavour disconnected piece Dgg. Both quark and antiquark must carry the same flavour. 



In what follows we will employ a Mellin transform in to undo the convolutions in eqn. ([T|): 
defining 

a{N,Q^)= [ dThT^~^a{Th,Q^), (3) 
Jo 

the factorization eqn. ([T]) becomes simply 

a(iV,Q2) = ^q(q2) ^ Dab{N;as{Q''))Fa{N,Q^)Fi,{N,Q^). (4) 

a,h=qi,qj,g 

Clearly Dab = D^a, and thus since the strong interaction is CP even Dg.g. = Dg.g., Dg^g^ = Dg^g^ 
and Dg^g = Dg.g. Thus if we define singlet and nonsinglet quark (plus antiquark) distributions and 
quark-antiquark luminosities as 

Q ^ J2^Q^ + Q^), (Q) ^ + (5) 

i i 

{QQ} ^ J^Q.Q* {QQ)^j^Y.''^Q'Q'^ 

i i 

where Qi = xqi etc, then we can decompose the factorization eqn. (j4]) as 

a = a0[Dgg-{QQ) + Dgg{Q)Q + Dgg{Q)G 

+nfDgg{QQ} + UfDggQQ + 2nfDggQG + nfDggGG] . (6) 

Inevitably this flavour decomposition is rather more complicated than in DIS (see ref. [11]) because 
here one must consider the flavours of three partons: the two incoming ones and the one from which 
the virtual photon is emitted. 

We now consider the behaviour of the perturbation series for the coefficients Dab{N; as{Q^)). At LO 
the only process is qq — > 7* and thus Dgg = 1 (fig(T^), while all the rest are zero. At NLO the channel 
Q9 ~^ 97* opens: Dgg = 0{as) (figE])- Processes with two gluons, two quarks, or two antiquarks in 
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(a) (b) 



Figure 2: NLO diagrams contributing to D, 




(a) (b) (c) 



Figure 3: NNLO diagrams contributing to Dgg at NNLO. 

the initial state can only occur at NNLO: Dgg = 0(0^) (figl3]), and Dgq = 0(0^) (figE^ andfSb). At 
this order there is also a "flavour disconnected" process contributing to qq ^ J*'- Dqq = 0{a1) (figllt). 
Flavour disconnected processes contribute also to qg — > X7* [Dgg = 0{a\)^ figSt) and qq — > X7* 
{Dqq = 0(af), figlSt): note that these each have a — > qq^y* subprocess. It is easy to see that all 
these flavour disconnected pieces are two powers of down compared to the corresponding flavour 
connected piece, and that they must be proportional to nj. 

Now consider the behaviour of the coefficient functions at high energy, that is as ^ 0. Only 
some of the higher order corrections contain high energy logarithms, and these are the ones we wish to 
isolate and compute. To see the how these logarithms arise, and thus deduce the general pattern, note 
that at LO the gg anomalous is singular, ^gg ~ while the qq anomalous dimension is not, 7^^ ~ a^. 
Thus the NLO corrections to Dqq due to gluon emission from a quark line (figdja) are 0{as)'- further 
gluon emissions from quark lines give contributions of 0{a^). Similarly the flavour disconnected piece 
Dqq has no high energy logarithms. In fact both the terms in eqn. ^ in which an incoming quark and 
antiquark annihilate to produce the vector boson are completely regular as A ^ 0, in the same way 
that the nonsinglet contribution to DIS structure functions is regular: 

Dqq = 1 + C7^(|^2 + o{al), Dqq = 0{al). (7) 
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(a) (b) (c) 



Figure 4: (a) and (b) are NNLO diagrams contributing to Dqg, which introduce high energy logarithms, 
while (c) is a N'^LO diagram contributing to the flavour disconnected piece Dgg. 




(a) (b) (c) 



Figure 5: (a) and (b) are NNLO diagrams contributing to Dqg, while (c) is an N^LO diagram con- 
tributing to the flavour disconnected piece Dgq. 



However the NNLO contribution to Dgg can be singular, due to gluon emission from the incoming 
gluon line (figSh- andlDa), and is thus (D{asj^)- Further gluon emissions from the incoming gluon 
give extra powers of j^, while gluon emissions from a quark line just give an extra Os, so the most 
singular contributions are Oiosijf)''), with A; > 0, i.e. NLLx. Similarly the N^LO contribution to 
Dgg is 0{a'^j^), due to gluon emission from one of the incoming gluons: in general the most singular 
contributions here are C'(a^(^)^), k > i.e. NNLLx. Finally consider Dgg-. at NNLO besides gluon 
emission from quark lines, one of the initial quarks can turn into a gluon, the rest of the process then 
proceeding as for Dgg (figl5^ and [6)3). Since 'jgg ~ Ujff' these graphs are also singular, 0{asjf). 
Further emissions from the gluon line then give contributions of 0{as{jf)^), k > 1, and are thus again 
NLLx. In fact these graphs are at high energy identical to the corresponding graphs for Dgg, save for 
the extra factor of ^ due to the incoming quark rather than incoming gluon, so there is a colour-charge 
relation [11] between the leading singularities of Dgg and those of Dgg. Similar considerations apply to 
the disconnected contributions Dgg and Dgg-. the most singular contributions to Dgg are 0(a^(^)^), 
k > 1, and related through a factor of to those of Dgg (see fig. Ht), while those of Dgq are 0(a^(^)'^), 
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k >2, and related through a factor of (^)^ to those of Dgg (see fig. Eb). 

Summarising, in the high energy limit the leading singularities of the coefficient functions have the 
schematic structure 

oo 
k=l 

oo 
k=l 

oo 

D,,iN,as) = {%LfY.^{al{^f), (8) 

k=2 

oo 
k=l 

oo 
k=l 

It follows that while the quark-quark and quark-gluon coefficient functions are NhLx, the gluon-gluon 
coefficient function and the disconnected contributions to quark-quark and quark-gluon are all NNLLx, 
and thus need not be considered in a NLLx calculation (though they may still be numerically large, 
since the gluon distribution is rather larger than the quark distribution). Our strategy in the rest 
of this paper will thus be to use /cy-factorization to compute the leading high energy singularities of 
the Dqg coefficient function to all order in perturbation theory, and then deduce those of Dgg using 
colour-charge relations. The high energy singularities are found by considering the off-shell process 
g*q 7*g calculated from the graphs in fig. [2] with the incoming gluon off-shell, in just the same way 
as one determines the coefficient function C2g in DIS using the off-shell process 7*5* qq [11]. The 
computation is complicated by the fact that just as in the DIS calculation there are collinear singularities 
as well as high energy singularities, so one has to take care to factorize the collinear singularities before 
extracting the high energy ones. In the next section we explain in some detail how this can be done. 

The hadroproduction of vector bosons, either or Z, has a very similar structure to Drell-Yan: 
in Mellin space the cross section can be written as 

a''{N,m^y) = a^ ^ D^,{N;asiQ^))Fa{N,Q^)F,{N,Q^) . (9) 

a,b=q,q,g 

where now 

aX = ^/^Gf (10) 

and Gp is the Fermi constant. The decomposition eqn. ([1]) is the same up to threshold effects, except 
that for Z-production ef in eqn. ([5]) is replaced by vf + af, where Vi, Oj are the vector and axial couplings 
of the Z bosons to the quarks, while for production 

{Q) = Q, (Q0) = J^l^^i.fQiOj, (11) 
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where Vij is the CKM matrix, Y^jVijV*j^ = 5ik, so the disconnected pieces of the quark-quark and 
quark-gluon coefficient functions may be simply added to the connected pieces. The perturbative QCD 
corrections to vector boson production are the same as the Drell-Yan ones at NLO, and they only begin 
to differ at NNLO because of the diagram with two incoming gluons and an internal quark loop, which 
gives an O(a^) contribution to the Z boson cross-section but vanishes in the case of a virtual photon or 
W^. The NLLx singularities are thus determined by the off-shell processes g*q ^ W^q' and g*q ^ Zq 
at 0{as) and are the same as those for Drell-Yan: the NNLLx singularities for Z production would 
however receive an extra contribution. 

3 High energy factorization 

The general structure of high energy singularities outlined in the previous section arises from the kx- 
factorization formalism [10], and in particular its extension by Catani and Hautmann to deal with 
situations in which the hard cross-section contains collinear singularities which must be factorized 
consistently [11]. In this section we explain how this works for Drell-Yan processes. 

We saw in the previous section how the singular contributions to the cross-section eqn. ([6]) are of 
essentially two kinds: those in which the electroweak boson is emitted directly from one of the incoming 
quarks or antiquarks, and those (the 'disconnected' pieces) where it is not. While the former can only 
acquire high energy logarithms through gluon emission from the second parton, the former acquire high 
energy logarithms through gluon emission from both partons. The dimensional regularised form of kx- 
factorization for the singular part of the Drell-Yan cross-section is thus (in dimensional regularization) 

as{N,Q') = ao(g2)||^^!-!!^S,,(iV,^;a,,e)(QW(iV,/.,e)) 

(.Ff {N, a., e)Q(°) (N, ^, e) + .^f {N, a„ e)^^ {N, ^, e)) 

J — 

(.Ff {N, |; a„ e)Q(°) {N, ^, e) + {N, a„ e)^^ (iV, e)) 

fi, e) + {N, ^; a^, e)G(°) (iV, /i, e)) I , 

(12) 

where 'Eqg{N,^;as,£) is the partonic cross-section for qg — > 7*5 with the incoming gluon off-shell 

k^ 1^.2 

(given at leading order by the diagrams fig. El and thus 0{as)), T,gg{N, ^;as,e) is the partonic 
cross-section for gg — > 'y*qq with both the incoming gluons off-shell (given at leading order by the 
diagrams fig. and thus 0{a'^)), J^f^ and !Ff^ are the bare quark and gluon Green's functions, and 
Q^^^ and G^^^ the bare integrated quark (plus antiquark) and gluon densities^ At LLx the quark 

^Note that our notation is slightly different to that in [If], since in particular we prefer to deal throughout with 
integrated parton densities. 
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Green's function is 



A'Hn, g;a„e) = gf f^f (iV, g;a„e) - nk' 6'~'-'^Hk)] , (13) 



while the gluon Green's function is 



The gluon-gluon transition function (containing the colhnear poles) is given by 

Tgg{N;as,e) =expl- —jgg{N,a)\, (15) 

where Se = exp[— e(— 7^; + ln47r)] and Ri'Jgg) is a scheme-dependent normahzation factor, regular as 
e — > 0: in the MS scheme 

m ^ f r(i - 7»)xo(t„) 1 / ,„> , T". '^■'(i)-i»'(i-7) 1 

''<^»»'=|r(i+.,,)[-7„xi(7„)il -i-h'*<'>^/„ — — r 

In the high energy limit 

7<,p(iV,a,)=74t)+0(«47f)') • (17) 



where Is {w) = + O j is the naive dual to the BFKL kernel xo- 

We can use the the /cr-factorization eqn. (jl2p to extract formulae for the bare coefficients in the 
expression eqn. Q at high energy in terms of the off-shell partonic cross-sections T,qg and S^^. For 
instance, comparing the coefficients of we have 



D(o)(iV,Q2;a,,e 



cf-^'k 
vrk^ 


^gg{N, 


k2 


(f-^'k 

^k2 


^gg{N, 





2\ 799 



= y ^^^9s(^'§7;as,e) ( jggR{7gg)Tgg{N;as,e), (18) 

where in the second line we substituted the gluon Green's function (jl4p . Introducing ^ = k^/Q^, and 
defining the impact factor 

V(7V,M;a,) = / ^^^^^"'^^^(iV,^; a„ 0) (19) 

JO 

we can write this as 

D(°)(iV,Q2;a„e) = (^^Y'' h,g{N,^gg;ash-g'Rhgg)Tgg{N-as,e) + 0{al{^)''). (20) 

Note that we can set e = in the impact factor, since collinear singularities are now regularised by 
M 7^ 0. The reason for the extra factor of in eqn. (|19p will become clear shortly. 
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Similarly, comparing coefficients of {Q^'^^)Q^^^ we find 

D^°^{N,Q^;as,e) = ^\ (^V'' h,g{N,-fgg;as)7ggR{7gg)rgg{N-as,e) 



Oil 

Ca 



DfhN, a,, e) - ^qg{N, 0; a„ e)! + 0(a2(^)'=), 



(211 



the second term (which is also singular) coming from the delta- function subtraction in eqn. 
Similar expressions for the bare coefficients Dqq, Dqg and Dgg may be extracted in terms of double 
Mellin transforms of S^^: we will ignore these in what follows as they are O(a^). 



4 CoUinear factorization at high energy 

We now need to factorize the collinear singularities in order to absorb them into renormalized parton 
distribution functions. In the case of a collinear safe process such as heavy flavour hadroproduction 
[14, 15] this is relatively straightforward: all the e poles are in the Green's functions, and in the MS 
factorization scheme are absorbed into the bare gluon distribution, i.e. G = TggG'^'^\ However when 
we consider a process in which the partonic cross-section is not collinear safe, such as DIS or DY, the 
factorization of the collinear singularities is rather more involved [11]. 

Firstly, we must separate out the valence and nonsinglet quark distributions from the singlet, since 
the former factorize multiplicatively while the latter mixes with the gluon. Defining 

Vi = Qi-Qi, Ni = Qi + Qi-^Q (22) 
so Ni = 0, we can rewrite the various quantities appearing in eqn. ([6]) 

i 

{QQ} = l^m + ^Qf-V^) = l{N'}-l{V'} + ^Q^ (23) 

i 

i 

When written in terms of N, V, Q and G the Drell-Yan cross-section eqn. ^ separates into three 
distinct pieces: a = a^^ + + a^^ , where 

CT^^ = ao[lDqg{{N')-{V')) + lnfDqq{{N'}-{V'})], (24) 
= ao[iDqq + ^Dqg){N)Q + Dqg{N)G], (25) 

(7^^ = ao[{Dqq + ^Dgq)QQ+DqgQG 

HnfDqq + ^Dqq)QQ + 2nfbqgQG + UfDggGG] . (26) 
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Now the factorization of collinear singularities for the valence and nonsinglet distributions Vi and 
Ni is multiplicative: 



Vi{N,as)=r^JN,as,e)Vl^\N,as,e), N^iN,as) = r^JN,as,e)N^^> iN,as,e) 



(27) 



so the nonsinglet cross-section a^'^ factorizes multiplicatively. Furthermore the transition functions 



and contain no high energy singularities: at high energy the both reduce to unity up to corrections 
of 0{aa)- However singlet quark Q and gluon G mix: 

Q = r,, + 2n fV^g G = g(o) + Tgg , (28) 

and the transition functions now contain high energy singularities. These were determined to NLLx 
in [11]: for example 

-^qg{N,a)Tgg{N,a,e), (29) 



Tqg{N, as.e) 



1 

e Jo 



a 



where Tgg is given by eqn. p^ . and thus Tqg contains (singular) terms of 0{as{j^)^). Tgg and Tqq are 
given in terms of Tgg and Tgg by colour-charge relations: 



r^, (iV, as,e) = ^ [Tgg (iV, a„ e) - 1] 



while writing Tgg = 1 -|- 2nfTg^ , 



rf/(7V,a„e) 



qq 



£ Jo 



aaSe 



da 



a 



hqq+7qqrgg] = g [r,,(iV,a„.) - ^^'ggiN)] 



(30) 



(31) 



since at high energy jgq{N,as) = ^[ygg{N,as) - as'jqgiO)], where jggiN^as) = as'jggiN) + 0{a1). 

Thus if we consider the collinear factorization of the singlet cross-section , substituting eqns. (|27l28p 
into eqn. ([2^ . the gluon and quark contributions to are given respectively by 







qg 



(Dqq + ^Dgg)2njTgg + DggTgg] {N ) G^''^ , 



^0 '^q 



so we have the collinear factorizations 



-.AT 
qq 



(Dqq + 2h^qq)rqq + DggTgg] {N^^^)Q 



(0) 



^(0) 

qg 

^ qq ^ 2nf ^qq > 



{Dgg + 2n^Dgg)2nfTgg + DggTgg 

{Dqq + i^Dgq)Tqq + DqgTgq 



(32) 
(33) 

(34) 
(35) 



Now consider first the gluonic contribution. Using eqns. ([7]) and eqns. ([8]), then since Tgg = 1+0(0^), 



while Tgg = 0{as{^f) 



Df^=Tgg+DggTgg+0{al{^f 



(36) 
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We would like to factorize out an overall factor of ^gg, in order to compare to the /cj^-factorization 
eqn. (jl2p . This can be achieved by taking the logarithmic derivative with respect to [11]: since 



d 



91nQ2 
d 



Tgg{N,as{^r,e) = ^gg{N,a,{^f Se)Tgg{N,as{^Y ,e) 



^^T,g{N,as{^y,e) = ^,g{N ,as{^f S,)T gg{N ^a^i^f ,e), 



we find 



d 



d\nQ 



2 99 



799 + 799^99 + ^« 



_d_ ' 



On the other hand, from eqn. (120p the high energy factorization gives 



d 



dlnQ 



2 99 



hgg{N, -Igg- U s) R{l g g)T g g{N ' U s , S) + 0{al{^f) 



(37) 
(38) 

(39) 
(40) 



Comparing eqns. (|39p and (|40p . we can factor all the collinear singularities into Tgg on each side, and 
thus we must have 



7,3(iV, as) + ^s{^)Dgg{N, a,) = h.giN, 7^; as)Rils) + 0(«'(t)') ' 



(41) 



where we also used eqn. (fT7|l . We can thus reduce the calculation of Dgg in MS factorization to a 
calculation of the impact factor eqn. ([19]) since 7^, 7^^ and i?(7s) are already known to the required 
order [11]. 

The quark coefficient function Dgg may be determined similarly: since Dgg = 1 + (D{as) and 



= 1 + 0{as), eqn. ([35]) may be written 



DiS = ^99 + rf/ + 1^,9 r99 + (^)') • 



(42) 



At high energy this is given by eqn. (|2ip. The subtraction term is just the on-shell partonic cross- 
section, i.e. the bare coefficient function evaluated at 0{as)- It may be evaluated using eqs. (|36p and 
([2H]): 

S,,(iV,0;a„e) = ^7°,(A^) + «sI)°<,(iV) + O(a^) (43) 

where the first term is the collinear singularity from the quark-gluon splitting and the second the 
renormalized coefficient function to 0{as)- Substituting eqs. ([36]) . ([12]) and ([^3]) into the bare colour 
charge relation eqn. (|2ip we obtain: 



Dgg + + Dgg Tgg = ^ [V gg + Dgg Tgg - - a.Z)",] + O {^f) . (44) 



Using the colour charge relations (I31|) and (|30p all the singular terms cancel as they should, and we are 
left with the renormalized colour charge relation for the Drell-Yan coefficient functions 



Dgg {N, as 



% {Dgg (N, as) - asD^iO)] + 0(a^ (ff )') . 



(45) 
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Thus the high energy behaviour of the pure singlet quark coefficient function Dqq may be determined 
directly from that of the quark gluon coefficient function Dqg found from the eqn. (jl9p through eqn. (|4ip . 

Finally we consider the factorization of the remaining contribution to the Drell-Yan cross-section, 
ijSS given by eqn. (j26p . The terms in the first line are 0{as{^)^), while those in the second are 
C'(a^(^)'^), so we consider them separately. Prom the coefficients of Q^'^^G^^^ and Q^'^-'Q^''-* in the first 
line we now find (in the same way that we derived eqns. (|34p and (|35p ) 

Df^ = 2{Dqq + ^Dqq)Tqq2nfTqg + Dqg{TqqTgg + 2nfTqgTgq), (46) 

Keeping only terms of 0{as{j^)^), these simplify to eqns. (f36l) and (|42]) . as they must. 

From the second line of eqn. ()26p consider the gluon-gluon term: it is easy to show that this 
may be used to determine the resummed coefficient function Dgg in terms of an impact factor com- 
puted as the double Mellin transform of the second derivative of the off-shell partonic cross-section 
gi, ^;as,e) in eqn. ([T2|): the result is 

jl+J,g7sDgg+^',Dgg = h g g{N , ^ s , 1 s'. ^ s) R{l sf + 0{al{^f) , (48) 

where 

hgg[N, Ml, M2; as) = MfMi / dCi / ^^2 C^''-%giN, ^1,^2; Os, 0). (49) 

Jo Jo 

However in this case the bare colour-charge relations implicit in the A:T-factorization do not lead directly 
to simple colour-charge relations among the coefficients Dqg , Dqg and Dgg , as they did for heavy quark 
production [14], because here there are also contributions of the same order from 0{a1{j^)'') terms in 
Dqq and Dqg, and these cannot at present be computed until /cT-factorization is extended to NNLLx. 

For the remainder of this paper, we will work entirely at NLLx, computing the impact factor hqg 
and thus all the C'(as(^)'^) contributions to the coefficient functions Dqq and Dqg. 

5 The off-shell cross-section 

In this section we compute the off-shell cross-section Sgg(r, a^, 0) for the process 

g*{k) q{p) 7*(g) q{p') , 

where 

k = Xipi+k, P = X2P2, 

q = zixipi + (1 - Z2)x2P2 + q, (50) 
p' = (1 - zi)xipi + Z2X2P2 k - q, 

and we define the dimensionless variables 

02 |lc|2 

r = ^, f=|L, ,51) 
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with = and v = 2j;iX2Pi • P2- We need to consider two diagrams, given in fig. [21 Since soft 
and collinear divergences are regulated by tlie off-sliellness of the incoming gluon, we may perform the 
whole calculation in four dimensions. The two-body phase-space is then given by 

^^^'^ = (2^)3 '^(g' - Q'np")i^^rs^'\k+p- q-p') 



6iq'-Q')S{p") 



(2^)2 

- ^dzidz2d\6{{l- Z2)ziu-\qf -Q^)6{{l-zi)z2iy-\k-q\'^). (52) 

The squared matrix element is computed using the usual eikonal polarisations for the gluon; the photon 
indices are contracted with g^'^ , because of conservation of electro-magnetic current, just as in the on- 
shell calculation [1]. The result is 

(Q2 + |k|2)(|k|2 - 2k • q) + 2(k • q)2 + \k\s - t)] } (53) 
where the Mandelstam invariants are 

s = (p + kf = u - m"^, 

t = {p - qf = - ziv, 

u = Q2_|k|2_s_i. (54) 

In the on-shell limit |k| 0, when averaged over the angle •& between k and q, the off-shell matrix 
element eqn. (|53l) reduces to 

i2\ _ ^Idl (t , s _ ^Q'^u- 



2 

^7t 



hm {\M\% = - + - + 2:^ , (55) 

|k|^0 Nc t St i 

in agreement with the standard calculation. 

To perform the phase space integration we first change the variable of integration by introducing 
A = q — zik, and we then use one delta function to perform the integral over ^2, and the second to 
perform the integral over = |qp — 22;ik • q+ zf|kp, i.e. we write 

1 dzidz2 ,2. .A-,,, i„i2 "'^ 



da>(2) = — / \ d^A 6{{l - Z2)ziu - |qp - Q^)6{Z2 



= j^dzidi)dz2dA^6{z2-^E^)S{A^-{l-z,)[zi{u-\k\^)-Q^]). (56) 

The remaining integral over ■!? then runs from zero to 2tt, while the integral over zi runs from Q'^ jiy — 
|k|2) to one. The dimensionless cross-section can then be written as 

n T l'^ /•27r J O 

S,,(r,e;a.,0) = ^-y_^dzi -\M\\,^^^,, (57) 



i-t5 
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where ^2 is the squared matrix element (|53|) with the two (5-functions (|56p imposed. The 

normaHzation is consistent with eqn. ([2]), with in particular the factor has been absorbed into the 
LO cross-section. The two remaining integrations can now be performed explicitly: the result is 



(l-T)(l-gr) 



ln( ^-'y^ )pi(r,g) 



+ 



(58) 



where the polynomials 
P2{r, 



r2 + (1 - r)2 + r^^ (10 + ^ + ISt^C - 6r(3 + 2^) , 
-1 + 36r^^3 + 7r(2 + ^ - 6r^C^(15 + 7^) 

+2t3^(35 + 49^ + 4^^) - t2(15 + 71^ + 14^2)_ 



(59) 



Note that this cross-section has been already averaged over the azimuthal angle of the incoming gluon, 
consistently with eqns. (jlSp and ()19p . where the Mellin integral is taken with respect to ^ = k^/Q^. 
The on-shell limit ^ ^ is 



Sgg(r,^;as,0) = — 



rP,,(r) In (i^)+lr(-i + 7r-fr2) 



(60) 



where Pqg{T) = Tr{t'^ + (1 ~ t)'^) is the LO g(7-splitting function. Because of the off-shell regularisation, 
this result is not the same as the usual MS NLO coefficient function [1]: 



rPqg{T) ( ± + ln47r - + In {^^) ) + \t{\ + 3t 



(61) 



though the collinear singularities match as they must, with ln(l/^) 1/e + ln47r — 7^;. 



6 The quark-gluon impact factor 

The quark-gluon impact factor is defined as the double Mellin transform of the double logarithmic 
derivative of the off-shell cross-section (see eqn. (|19p : 



l^co pi 

hggiN,M;as) = M^ ^'-^ T^-iS(T,e;a„ 
Jo Jo 



0) 



(62) 



Because the off-shell cross-section eqn. (I58p contains a factor Q{t ^ — ^ — 1), it is useful to change the 
integration variables according to 



a = T^, (3 



(63) 
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the Jacobian determinant is and the 0-function condition is satisfied for all a, /3 G [0, 1], so 



271" Jo Jo 



In 



l-(l-a)/3 



1 (1 - /?) 



(64) 



where the polynomials di and d2 are now 



di{a,P) = l + a^ -2[3 + l2af3-22a^l3 + l2a^f3 + 2f3'^ -22a(5^ 

+560^/3^ - 54a3/32 + ISa^/J^ , 
d2 (a, /3) = 1 - 6a + Sa^ - 14/3 + 71a/3 - QSa^/? + 42a^/? + 15/3^ 

-850/3^ + leOa^/?^ - 126a^/32 + 36a'' /J^ . 



(65) 



The integrals containing the logarithm can be further simplified through integrating by parts with 
respect to the variable j3: 



' da d/3 a^'-^+P{l - a)^-^^/3^-^^+'? In ^-^^ 



JO 

1 



l + N -M + q 



JO 



da d/3a^-i+?'(l-a)^-*^/3^-^+« 



1 - (1 - a)(3' 



(66) 



where p and q are integers, the exponents of a and (3 in each terms of the functions di{a, j3). Thus the 
only nontrivial integral which we need is 



da d/?a*^-i+P(l-a)^"*^/3^-*^+'? 
Jo 



1 



1 - (1 - a)/3 

oo „i „i 

da d(3a^'-^+P+''{l - a)^~M ^N^M+g+k 

^ r{l + N - M + k)r{M + p) 
^^T{1 + N + k + p){l + N - M + q + k) 

T{M +p)T{l + N - M) 
T{l + N+p){l + N-M + q) 
sF2{l,l + N - M,l + N - M + q;l + N + p,2 + N - M + q;l) 



(67) 



T\M+p) 



r(i + - M)r(i + - M + 



T{l + N + p)T{l + N + p + q) 
3F2{M +p,N + p,l + N - M + q;l+N +p,l + N +p + q;l), 



(68) 



where 3^2(0; b, c; d, e; z) is a generalised hyper geometric function, and in the last line we used Thomae's 
theorem. The final result for hqg eqn. (j64p is then a fairly complicated sum of terms involving the 
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generalised hyper geometric functions 3F2: 



hgg{N,M;as,0) 



as VV,.2t.2.,. I T{l + N-Mni + N-M + q) 



-p + q) 



p=0 q=0 

3F2{M +p,N +p,l + N - M + q;l + N + p,l + N + p + q;l) 



[1 + N -M + q) 



+ ^l+p,l+g 



(69) 



where 



( 1 


-2 


2 


\ 




/ 


1 


-15 


29 


-15 


\ 





12 


-22 





' 2 




-6 


77 


-156 


85 




1 


-22 


56 







8 


-106 


258 


-160 







12 


-54 










42 


-168 


126 




V 





18 


) 




V 








36 


-36 


/ 



(70) 



Note that in the colhnear hmit M — > eqn. ()69p is regular. In fact since only terms with p = survive, 
and 3i^2(0, N,l + N + q;l + N,l + N + q;l) = 1, 



a. 



hgg{N,0;as,0) = ^Tr 



+ 3iv + 4 



27r ""{N + 1){N + 2){N + ?,) 



(71) 



as it must be from eqn. (j4ip . This is a highly nontrivial check on eqn. (j69p . 

In the limit N ^ Q (i.e. in the high energy limit), the sum (j67p may be performed more easily, since 
the generalised hypergeometric function reduces to a sum of rational functions of F-functions: 



Using this in eqn. 



M n-M+q 



da / df3a^^-^+'P{l-a)-^^ 
Jo 

r{i- M + k)r{M + p) 

~ ^T{l + k+p){l- M + q + k) 



1 - (1 - a)P 



T{M + p) 



r(M + p)T{l -M - p)T{l -M -p + 



r(i + (?) 



p-1 



i=o 



Tjl-M-p + j) 
{1- M + q-p + j)j\ 



the result remarkably simplifies to 



' ' ' ^ 27r ^ (1-M(2-M) 3-M 



(72) 



(73) 



The same result may be obtained by setting = in the generalised hypergeometric functions in 
eqn. ([69]) . and simplifying the result. The infrared singularity at M = 1 is a triple pole, as expected 
from the counting of soft and collinear singularities [12]. 
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7 Drell-Yan coefficient functions at high energy 



The Taylor expansion of eqn. (j73p about M = gives 



hqgiO,M;as,0) 



—Tr- 
2tt 3 



11 



/85 



1 + —M+ — + 



6 



V36 



575 Il7r2 

^ M-" 

216 18 



1296 108 15 ' ^ ^ 



(74) 



In order to compute the high energy behaviour of the MS coefficient function, we need the scheme 
dependent factor: 



R{M) = 1 + f CsM^ - f C4M^ + O (M^) 



The BFKL anomalous dimension is 



Is 



ttN 



+ 2C3 



ttN 



+ 2C5 



ttN 



+ ... 



(75) 



(76) 



while quark anomalous dimension in MS may be found in [11]: 

2 



7,9 (A^, 



a. 



1 + 14 /^C'yl 

3 TT 9 



/82 



/122 25 
V243 ^ ¥ 



CaoA^ /146 14^ 



vr N 



+ 



(77) 



Substituting eqn. ([76]) into eqns. ([7^) and ([75]). and then substituting the results along with eqn. ([77|) 
into eqn. ()4ip . the leading terms on either side cancel, and the remaining terms give the NLLx expansion 



Dqg{N, as 



18^^^ 



29 



1+ — +27r^ — ^ 



6 



Ca Os 



+ 



/1069 11 



108 



vr 

CAas\ 



vr 



iv 



/9031 , 85 2 7 4 73 



648 18 



20 



vr N 



+ 



(78) 



The coefficients of O (a^) and O (a^) are in agreement with the high energy limit of the fixed order 
NLO [1] and NNLO [3,4] computations. This is a very non-trivial check of the procedure. The O (af) 
and subsequent terms are all new results. 

The high energy singularities of the quark-quark coefficient function are now easily deduced using 
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the colour-charge relation eqn. (|45p : we find 



Dqq{N, as) 





vr N 



) 



2 



+ ... 



(79) 



Again this result checks at O (a^) against the NNLO result of ref. [3,4], while the terms from O (q^) 
onwards are new results. 

8 Summary 

We have analysed the Drell-Yan (and thus also vector boson production) process in the limit of high 
partonic centre-of-mass energy, along the lines of the analysis of DIS presented in [11]. We showed how 
to factorise simultaneously high energy and collinear singularities from the hard coefficient function, 
computed the quark-gluon impact factor, and thus deduced the leading high energy singularities in the 
Drell-Yan coefficient functions in MS scheme to arbitrarily high orders in ag. Our results agree with 
the known results at NLO and NNLO, while providing new results at N^LO and beyond. 

It will now be possible to use these results to perform an all order resummation of high energy 
logarithms in Drell-Yan processes, using the techniques developed in [12] which have already been 
applied to deep inelastic processes [13]. A global analysis of the effect of high energy resummation 
on parton distribution functions and LHC benchmarks is thus now a definite possibility for the near 
future. 
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